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Abstract 

Necessary conditions of optimality are derived for trajectories whose 
structure is limited to climb, steady cruise, and descent segments. The 
performance function consists of the sum of fuel and time costs, referred to 
as direct operating cost (DOC). The state variable is range to go and the 
independent variable Is energ". In this formulation a cruise segment always 
occurs at the optimum cruise energy for sufficiently large range. At short 
ranges (500 n. mi. and less) a cruise segment may also occur below the 
optimum cruise energy. The existence of such a cruise segment depends 
primarily on the fuel flow vs thrust characteristics and on thrust constraints. 
If thrust is a free control variable along with airspeed, it is shown that 
such cruise segments will not generally occur. If thrust is constrained to 
some maximum value in climb and to some minimum in descent, such cruise 
segments generally will occur. Computer calculations of typical short-range 
trajectories obtained about a 1% cost penalty for contraining the thrust. 

Index categories: Flight Operations; Guidance and Control; Navigation; 

Communication; Traffic Control. 
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Nomenclature 
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P 

S FC 
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T ,T, 
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* fuel coat factor, dollars/pound (dollars/kilogram) 

* time cost factor, dollars/hour 

- drag force 

* first and second partial derivatives of drag with respect 
to airspeed 

- cruise distance 

* desired distance to fly 

■ total climb and descent distances, respectively 

* total aircraft energy in units of altitude 

* cruise or maximum energy 

- optimum cruise energy 

* initial and final energy 
« rate of change of energy 

* acceleration of gravity 

* Hamiltonian, dollars per unit of energy 

* altitude in feet (meters) 

a components of the Hamiltonian 

= value of performance function in dollars or pounds (kilograms) 

* operands under the minimization operator in H 

* lift force 

* integrand of cost function or cost per unit time 
= thrust specific fuel consumption per hour 

- n c h partial derivatives of S^ c with respect to (•) 

* thrust in pounds (kilograms) 

" climb and descent thrusts, respectively 


2 



t 


time 


t 


c 


t 

t 


d 

f 


V 

V 

c 

v .v, 

up’ dn 

V 
w 

V V 
vup* wdn 


U 

W 

w 


f 

f 


x 

x *Xj 
up dn 

3 


Y 

AR 

AT, AV 
A(E c ) 


7T ,H 

up dn 

<P 

A 


A 

opt 


■ time at end of climb 

■ time at start of descent 

■ total mission time 

* true airspeed 

* cruise speed 

■ climb and descent airspeeds 

* wind speed along flight path 

■ wind speeds in climb and descent segments, respectively, 
functions of altitude 

* aircraft weight in pounds (kilograms) 

* total mission fuel in pounds (kilograms) 

■ fuel flow rate in pounds per hour (kilograms/hour) 

■ distance flown in nautical miles 

- climb and descent distances, running variables 

■ parameter defining direction of control perturbations 

■ flight-path angle, radians 

“ length of control perturbation 

* thrust and speed perturbations relative to cruise conditions 

■ cruise cost at cruise energy E^, dollars per nautical mile 

* throttle setting 

* throttle settings in climb and descent, respectively 

- costate variable 

= cruise cost per unit distance 

* optimum cruise cost over all energies, per unit distance 
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I. Introduction 


RECENT research on aircraft trajectory optimization has focused on the 
structure of optimum cruise trajectories. Steady-state cruise is generally 
hot optimum for minimum fuel performance (Ref. 1); but the performance 
improvement relative to steady state cruise is unknown, because the actual 
nonsteady optimum cruise has not been computed. However, if the steady-state 
cruise satisfies first-order necessary conditions, Speyer (Ref. 1) shows, 
in an example, that the performance improvement of a particular (though 
nonoptimum) cyclic cruise is about 0.1%. This improvement, if representative 
of the optimum cyclic cruise, is not economically significant. Nevertheless, 
the determination of the optimum cyclic cruise poses an interesting and 
unsolved problem. 

Even if economically significant, cyclic cruise could not be used in 
airline operation because it is incompatible with existing air traffic 
control procedures, disconcerts passengers, and decreases engine life because 
of increased cycling. Optimum trajectories, to be compatible with typical 
airline practice, should consist of: a climbout, a steady-state cruise, and 

a descent. Thus, at least for commercial airline applications, the optimum 
trajectory must be selected from a set of trajectories that is limited 
a priori to such types. 

A formulation of the trajectory optimization problem that constrains 
the admissible trajectories to those containing climb, steady cruise, and 
descent was given in Refs. 2 and 3. In this formulation, energy height was 
used as the independent or timelike variable in climb and descent, thus 
forcing energy to change monotonically in these segments. This formulation 
of the constrained trajectory optimization problem is also adopted in this 
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paper. The use of energy eliminates the Integration of a separate adjoint 
differential equation and simplifies the calculus of variations problem to 
one requiring only pointvise extremization at algebraic functions (Refs. 2 
and 3). 

While the problem formulations here follow that of Refs. 2 and 3, this 
paper extends those earlier results and places them on a firm analytical 
foundation through the use of the calculus of variations. Also, an evaluation 
of the constrained optimum trajectories by airline operators indicated an 
interest in the additional constraint of setting the thrust to some maximum 
during climb and to idle during descent. An examination of this procedure 
raised the following questions. How do the constraints on thrust and, more 
generally, the aerodynamic and propulsion characteristics affect the structure 
of the trajectories? Under what condition is the constrained thrust procedure 
optimum? What performance penalty is incurred by the constraint on thrust? 

The investigation of these questions begins in Sections II and III with 
the optimal control formulation of this problem and the derivation of the 
first order necessary conditions. Sections IV, V, and VI contain the details 
of the analysis relating the characteristics of the trajectories and the 
optimum controls to the propulsion model. Section VII contains a brief 
discussion of computed optimum trajectories for the Boeing 727-200 equipped 
with JT8D-7 turbofans. 


II. Optimal Control Formulation 

As stated in the introduction, we assume at the outset that the optimum 
trajectories have the structure shown in Fig. 1. This structure consists of 
climb, cruise, and descent segments, with the aircraft energy increasing 
monotonically in climb and decreasing monotonically in descent. Neglecting 


flight-path angle dynamics and weight loss due to fuel burn, the point mass 
equations of motions for flight in the vertical plane are 

(l/g)(dV/dt) * [(T - D)/W] - sin y (1) 

dh/dt ■ V sin y (2) 

dx/dt ■ V cos y + V — V + V (3) 

' w w 


with the contraint L ■ W cos y. In airplanes, unlike rockets, the rate of 
change of weight loss due to fuel burn introduces negligible dynamic effects 
in the trajectory optimization. Nevertheless, the effect of weight loss on 
a trajectory is important but can be accounted for without adding another 
state v Triable by continuously updating the weight as the trajectory is 
generated. If energy is defined as 

E = h + (l/2g)V 2 (4) 

then the familiar relation for the rate of change of energy is obtained by 
differentiating Eq. (4) with respect to time and substituting the right-hand 
sides of Eqs. (1) and (2) in place of dV/dt and dh/dt, respectively: 

E H dE/dt = [ (T - D)V]/W (5) 

The cost function to be minimized is chosen as the direct operating 
cost of the mission and consists of the sum of the fuel cost and the time cost: 


J - 


c,W, + c t 
f f t t 


( 6 ) 


where c^ and c^ are the unit cost of fuel and time, respectively. Setting 
c^ = 0 results in the familiar minimum fuel performance function. In integral 
form, the cost function becomes 
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P dt 


(7) 


J-/' f CVf + 't> dt = /' 

0 0 

It 1s assumed that the time to fly t^ is a free variable, but the 
distance to fly is a specified quantity d^. Following the formulation in 
Ref. 2, we now write the total mission cost as the sum of the costs for the 
three segments of the assumed trajectory, illustrated in Fig. 1. 


* 

/ 


P dt + 


(d. - d - d, )A 
f up dn 




P dt 


( 8 ) 


climb cost 


cruise cost 


descent cost 


where X designates the cost at cruise energy E^. Next, we transform the 
integral cost terms in Eq. (8) by changing the Independent variable from time 
to energy, using the transformation dt = dE/E, 


/' 


< p/E W dE + (d f - d up - d dn n + 


/ 

E, 


(P/|E|]^ 0 )dE 


(9) 


The transformation uses the assumption that the energy changes monotonically 

* 

in the climb and descent. This places strict inequality constraints on E, 
as shown in Eq. (9). Also in Eq. (9), the integration limits have been 
reversed in the descent cost term. In this formulation the cost function is 
of mixed form, containing two integral cost terms and a terminal cost term 
contributed by the cruise segment. 

With the change in Independent variable from time to energy, the state 
equation (Eq. (5)] is eliminated, leaving Eq. (3) as the only state equation. 
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Furthermore, we note that the performance function [Eq. (9)} depends on the 

distance state x only through the sum of the climb and descent distances 

d + d, . Therefore, the state equation for the distance Is rewritten In 
up dn ’ 

terms of this sum: 

d( *up + X dn ,/dE ■ (v up + + < v dn + V wdo )/ l E l ) E<0 < 10 > 


Here, the transformation dt ■ dE/E was used again. Also, Eq. (10) provides 

and 


for Independence In the choice of climb and descent speeds V and V 


dn 


th.r wind velocities V and V , . Wind velocities in climb and descent 

wup wan 

are allowed to be independent of each other; generally, different wind 
conditions will prevail in physically different locations of climb and 
descent. The wind velocities can also be altitude dependent. The effect of 
altitude dependent winds on the optimum trajectories is discussed in Ref. 3. 

Necessary conditions for the minimization of Eq. (9), subject to the 
state equation [Eq. (10)] are obtained by application of optimum control 
theory. (See, e.g. , Ref. 4, page 71.) Then, the following relations are 
obtained for the Hamiltonian and costate equations, respectively: 


H = min 


V ,V. 
up dn 

TT , 7T 

up’ dn 


' 

V + V V, + V . 

< 

(*) + + * 

up wup | dn wdn 


K>0 \I^E<0 

> 

^E>0 ^E<0 

i 


(ID 


dip/dE - -[ 3H/3(x + x, )] - 0 

up dn 


( 12 ) 


The right-hand side of the Hamiltonian equation is minimized with 

respect to two pairs of control variables, one pair applicable to climb 

(V and 7t ) the other pair to descent (Vj and it, ). Since each term under 
up up dn dn 
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the minimization operator in Eq. (11) contains only one of the two pairs of 
control variables, the minimization simplifies into two independent 
minimizations, one involving climb controls, the other, descent controls. 

Also, since the right-hand side of the costate equation (Eq. (12)] is zero, 

\f> is constant. 

III. Transversality Conditions 

The transversality conditions are additional, necessary conditions that 

depend on the end point contraints of state variables (Ref. 4). The basic 

constraint in this problem is the range of the trajectory d^. However, 

d^ is a parameter in the transformed cost function, Eq. (9), and not a state 

variable. The final value of the state variable d + d , is, in this 

up dn 

formulation, subject only to the inequality constraint d yp + d^ n £ d f . This 

constraint is, of course, necessary for a physically meaningful result. 

This inequality constraint can be handled by solving two optimization 

problems, one completely free (d + d^^ < d^), the other constrained 

(d + d, = d,) and then choosing the trajectory with the lowest cost, 
up dn f 

Physically, the comparison is between a trajectory with a cruise segment, 

and one without a cruise segment. Considering first the free terminal state 

case d + d , < d,, we obtain the following relation for the final value 

up dn f 

of the costate ij/j 

<KE ) - 0(d, - d - d. )A]/(3(d n + d. )] - -A (13) 

c f up dn up dn 

This is the transversality condition for the free final state problem with 
terminal cost (Ref. 4). It shows that the constant costate value is the 
negative of the cruise cost. 
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Next, consider the case of no cruise segment. Then, the middle term of 
Eq. (9) drops out and the performance function contains only the integral corf- 
terms. This is the case of the specified final state d^ - d up + d dn J the 
corresponding transversallty condition yields <f»(E c ) “ In practice it is 

not necessary to compute the constrained terminal state trajectory if a valid 
free terminal state trajectory exists, since the addition of a terminal con- 
straint can only increase the cost of the trajectory. Therefore, this case is 
not considered further in this paper. 

In both cases the choice of costate determines a particular range. Since 
the functional relationship between these variables cannot be determined in 
closed form it is necessary to Iterate on the costate value in order to achieve 
a specified range d^. 

The last necessary condition applicable to this formulation is obtained by 
making use of the fact that the final value of the timelike independent variable 
E is free. Its final value is the upper limit of integration E £ in Eq. (9). 
Application of results in Ref. 4 provides the following condition: 

(« + tn« f - % - d dn )*(E)l/3E)) ■ 0 (14) 

which, when evaluated and sinq>lified, becomes 

{H + Id c <dA/dE) n E _ E " 0 (15) 

c 

where d is the cruise distance, 
c 

Condition (15) has the following physical interpretation. The value of 

the Hamiltonian H evaluated at cruise energy E fi is [after substituting 

Eq. (13) into Eq. (11)] the minimum increment in climb plus descent cost 

relative to cruise cost to make a ur.it increment in cruise energy. The 

product d (dX/dE) _ _ is the minimum increment in cruise cost for a unit 
C E"l 

c 

change in crufse energy. Condition (15) requires the optimum trajectory to 
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be such that the sun of these two increments be zero for a given cruise 
distance d £ and cruise energy E c . 

IV. Dependence of Optimum Trajectories on Range 

Equation (7.5) t together with knowledge of the salient characteristics 
of the crulBe cost A and the Hamiltonian H, can be used to determine the 
structural dependence of the optimum trajectories on range. 

Cruise cost at a cruise energy E £ and cruise speed is computed 
from the relation 

A(E c ,V c ) - [P(T,E c ,V c )]/(V c + V w ) ; constraints: |£ “ jjj (16) 

.'here the denominator is the ground speed in the flight-path direction. 
Examination of the term containing A in the relation for the performance 
function (9) shows that the value for A should be as small as possible at 
each cruise energy in order to minimize the total cost J. Therefore, the 
cruise speed dependence of A is eliminated by minimizing the right side 
of Eq. (16) with respect to V £ 

A (E ) - min P(T,E ,V )/ (V + V ) (17) 

C y c c c w 

c 

In this paper, A and are always assumed to be the optimum cruise 

cost and cruise speed, respectively, at a particular cruise energy E^. 

Except in high wind shear, the cruise cost as a function of cruise energy 
exhibits the roughly , arabolic shape shown in Fig. 2. For subsonic transport 
aircraft, the minimum of the cruise cost with respect to energy occurs close 
to the maximum energy boundary. This characteristic of the cruise cost 
prevails for essentially all values of the performance function parameters 
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and c fe . The quantities defining the optimum cruise conditions are E CO p C 
and A In Eq. (15), the derivative of the cruise coat function multiplies 
the cruise distance. Except under extreme wind shear conditions, the deriva- 
tive is mono tonic and crosses the sero axis at I ■ E 

c copt 

By distributing the minimization operator in Eq. (11) and substituting 
Eq. (13) in Eq. (11) , H can be decomposed into clisd> and descent components 
as follows % 


Hti.*(E c )l - 1 ♦ I dn 


(18) 


where 


up 


min 

V 

up 

ft 

up 


P - A(V + V ) 
up wup 


E] 


E>0 


I , ■ rain 

dn 

dn 

IT. 

dn 


P - A (V + V . ) 
dn wdn 


'E<0 


(19) 


In the preceding section, the Hamiltonian, evaluated at E = E^, was 

interpreted as the cost penalty to achieve a unit increase in errise energy. 

Extensive numerical studies of Eq. (18) for several comprehensive models of 

subsonic turbofan aircraft show H[E ,A(E )] > 0, for E £ E . Moreover, 

c c c copt 

the minimum cost penalty for increasing energy I is always positive and 
that for decreasing energy 1^ is negative; but the sum has never been 
found negative for models of currently used turbofans. While these character- 
istics have been established for several aircraft models, they are not 
intended to imply a generalization to all aircraft, since no physical laws 
prevent H from being negative. 

Consider first the case where H[E »A(E ) > 0. Then, Eq. (15) can be 

c c 


solved for the cruise distance d 
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d - -H(E « X(E )/ (dX/dE)_ « 
C C c E"E 


( 20 ) 


Since dX/dE < 0, but approaches zero as E approaches E the 

c copt 

cruise distance oust increase without limit as E approaches E ... Our 

c copt 

numerical studies have shown that the value of H tends to decrease as E 

c 

increases, but not enough to change this trend. Figure 3 shows the resulting 

family of trajectories, assuming H > 0 for all values of E • In this case, 

c 

Interestingly, nonzero cruise segments occur at short ranges and at energies 
below the optimum cruise energy E CO p t * Optimum cruise is approached 
asymptotically at long range. 

Consider next the case where H[E c »X(E c )] * 0. Then d £ ■ 0, i.e., 
no r vise segment is present for dX/dE < 0. However, Eq. (15) shows that 
d^ be nonzero if dX/dE < 0. This implies that for H * 0, cruise flight 

is optimum only at the optimum cruise energy E copt * Figure 4 shows the 
family of trajectories for this case. 

It is conceivable that both cases appear in the same problem, although 
this has not occurred in the examples we have studied. 


V. Thrust Optimization for Minimum Fuel Trajectories 
Evaluation of the Hamiltonian equation would be simplified if one of 
the two pairs of control variables, airspeed or thrust, could somehow be 
eliminated a p riori from the minimization. Since the pair of throttle 
settings, n and Ti^ n > is thought to be near its limit, we shall look for 
conditions where extreme settings of the throttle are optimum. The remainder 
of the paper examines only the minimum fuel case c^ « 1 and c^ ■ 0, with 
winds set to zero in order to simplify the derivation. However, the results 
can be extended to the more general cost function. 
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Por minimum fuel performance, che two terms in the Hamilton!*. » Eq. (19) 


1 n ■ min K , I , * min K 

U P it ,v up dB ^ ,v , dn 

up* Up dn* dn 


(21a) 


where 


up 


V f - XV 


["i.:. > 

(T- D)V vif 
up 

am w 

• K dn - 

[t - D|V dtt /W 


(21b) 


T< ’d n > <D 


An accurate nodal fop thrust and fuel fl°w generally Includes the 
functional dependencies! T(w,V,h) and &£(ir,V»h), In addition, these functions 
must be corrected for nonstandard temperatures and bleed losses. 

In previous work cm aircraft trajectory optimization (Ref. 5), a 
simpler model for fuel flow and thrust was used 

Wf - TS FC (V.h) i T Btn (V,h) S T < T^tt.h) (22) 


The critical assumption in Eq. (22) is independence of the specific fuel 
consumption Sp c from thrust. The virtue of this model lies in the insight 
it yields into the minimum fuel problem. If Eq. (22) is substituted into 
Eqs, (21b), one obtains 


S W 

rc 


up v, 


up 


T - (X/S Br )V 


T - D 
up 


T >D 
up 


dn 


hdt 

V dn 


T dn ~ (>/8 FC )V d« 


l T dn * D 


^dn <D 


(23) 
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For any fixed values of V or V* , the operand functions for the 

Up cm 

minimization and are hyperbolae with poles at T * D, The 

numerator aero mist be to the left of the pole on the thrust axis for energies 

less than cruise energy. Figure 5 is a typical plot of these functions. 

Clearly, maximum thrust minimises K and idle thrust minimizes K. for 

up dn 

any E < E c » proving that the limiting values of thrust are optimum for this 
propulsion model throughout the climb and descent trajectories. This result 
also implies that the departure from the extreme thrust values found for the 
more general propulsion model is directly attributable to the nonlinear 
dependence of fuel flow on thrust. Conversely, the need for throttle setting 
optimization can be determined a priori from the fuel flow vs thrust 
dependence for a particular engine. Such data are found in the engine 
manufacturer's performance handbook. 

VI. Evaluation of Hamiltonian at Cruise 

We have seen in a preceding section that the value of the Hamiltonian 

computed at cruise energy E c determines the structure of the trajectories 

near cruise. Here we shall relate the existence of cruise below E _ to 

copt 

specific engine and aerodynamic model parameters. This is done by substi- 
tuting truncated Taylor series expansions of fuel flow and drag as functions 
of airspeed and thrust into the expression for the Hamiltonian. The location 
of the minimum with respect to the controls as well as the value of H can 
then be determined as functions of the Taylor series coefficients at E * E^. 

How should one pick the point in the control space about which to make 
the expansion? Computational experience in Refs. 2 and 3 has shown that 
the minimum is in the neighborhood of the optimum cruise speed and throttle 
setting, corresponding to the given cruise energy. This suggests that the 
cruise controls should be picked for the expansion point. 


15 



The fuel flow end dreg functions expended to second order ebout the 


cruise controls T « T . V ■ V ares 

c c 



Vic + (Vrc T + s Fc) iT + T c s rc v AV + < l/2 H 2s »c T 

+ (1/2)(t c S TCtv + S FCv y V At + «/2)T c S AV 2 


+ T c S FC , ' 4T 
T 


) 


+ higher-order terms (24) 

2 

D - D(V ,E ) + D AV + (1/2)D „ AV + higher-order terms (25) 

C C V 

The subscripts to S_ r end D designate the partial derivatives with respect 
to the subscripted variable. Note that the expansion allows for a general 
fuel flow model in which specific fuel consumption can be thrust dependent. 

Before substituting Eqs. (24) and (25) into the expression for H, we 
observe that H is singular at cruise with T * and V = V^, because 
both numerator and denominator are identically zero at that point. Figure 6 
plots the loci of the numerator and denominator zeros of K up and K dn in 
the control space at E ■ E . It is proved in the appendix that the locus 
of numerator zeros is tangent to the locus of denominator zeros at the 
optimum cruise controls. For E < E £ the two loci have no points in common. 
The two loci can be tangent but cannot cross, since otherwise controls 
would exist that would make the Hamiltonian infinitely negative, a result 
ruled out as physically meaningless. 

Upon substituting Eqs. (24) and (25) into Eqs. (21) using the tangency 
condition (A4) derived in the appendix, the following expressions for K^ p 
and K. at cruise energy are obtained: 
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Terms above second order have been neglected since we are investigating a small 
neighborhood of the cruise point. Expression (26) represents K up if the quan- 
tity under the absolute value sign is positive and Kj n if it is negative. 

Since the cruise point at AT ■ 0 and AV * 0 gives the undefined value of 
0/0 for Eq. (26), it is necessary to evaluate the limit as AT and AV approach 
zero. If the limit exists it must be independent of the direction from which the 
cruise point is approached. lo compute the limit and investigate the neighbor- 
hood of the cruise point, a polar coordinate system centered at the cruise point 
is used to define control perturbations. Let AR and 0 define control 
perturbations AT and AV as follows: 

AT - (D v + 8)AV (27) 

AV » AR//l + (6 + D )" 2 , AT - AR(8 + D y )// 1 + (6 + D y ) 2 (28) 

The parameter 0 defines a direction relative to the reference direction of the 
line AT * AV. The reference direction 8-0 is excluded from the control 

space since it is along the direction of the locus of T ■ D at the cruise point. 
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After substituting Eqs. (28) into 8q* (26) and taking the limit of the 
resulting expressions as AR •* 0, one obtains for any MO 


limit 


(W/V C )(S FC 


+ Vkl 


J limlt 


- <-w/v c )< Sfc 


The limit is thus well defined, since it is independent of the approach 
direction in each region. ttowevet, it remains to be shown that the limit 
value is in fact the minimum of Eq. (26) with respect to the perturbation 
controls. This question is investigated for two cases, one where is 

independent, and the other where it is dependent on thrust. 

Case (A): S Independent of Thrust 

FC 

Along the direction defined by AV ■ 0, i.e., along the thrust direction, 

Eq. (23) can be used directly to determine the dependence of the functions 

on T and T, under the minimization operator. Since at V - V , 
up dn c 

D(V ,E ) « T - (A.S__)V , Eq. (23) reduces to 

C C C rl C 

K „p ■ (H/V c )S FC • K dn ' ( - M/V c )S FC (30) 

showing that at the cruise speed these functions are independent of 

thrust. This result is not restricted to small perturbations relative to 
the cruise thrust. Along other directions, the truncated Taylor series form 
[Eq. (26)] must be used. After setting to zero all thrust dependent 
derivatives and substituting Eqs. (28) into Eq. (26), the following expression 
is obtained: 
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( 31 ) 


*“«’ } hs fc 

±1 + 

K.‘l" + V + VrcJ 

|AR 

2 1 a | S Fr , /l + (0 + D v ) Z 

w ‘ - (V + AV) 

K . C 

dn 


D 2 AV 

1 V — j 


2 1 3 1 /i + (8 + D v ) Z 


where the positive sign applies to and the negative sign to Kj ft . 

The characteristics of these functions depend on the drag and specific fuel 

consumption derivatives. The drag derivatives D and D 2 are both 

positive, since the aircraft will certainly operate on the "front" side of 

the thrust-required curve. The dependence of S pc on speed for a typical, 

currently InBervice turbofan engine at cruise energies exhibits a slight 

upward curvature above Mach 0.4, as shown In Fig. 7, implying that both 

S_ and S__ are positive in the range of interest between Mach 0.4 
FC V V2 

and 0.9. The slight curvature of S„ r indicates that a quadratic function 

r L 

can accurately model the Mach number dependence of in the Mach range of 

interest and not just in a small neighborhood of the expansion point. Also, 

at typical cruise conditions one finds that D 2 > (2Sp C + T c Sp C ). 

v V y2 

Therefore, for any 0, the denominator of Eq. (31) goes to zero before the 
numerator does as AR is increased from an initial value of zero. Moreover, 
the slope of the operand function with respect to AR increases as $ 
approaches zero. The effect of AV can be neglected, since >> AV. 

These observations lead to the conclusion that the functions in Eq. (31) 
slope upward in all directions as AR increases, except In the direction 
parallel to the thrust axis, along which the slope is level. Figure 8 shows 
a family of plots of the operand functions as 8 varies over its range. The 
limiting values of these functions at the cruise point (±W/V c )Sp c are 
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therefor* also the global minimum*, and the value of the Hamiltonian, which 
li the sum of the two components, is zero. At the cruise energy, furthermore, 
the optimum climb and descent .speeds are equal to the optimum cruise speed. 

The optimum climb and descent thrusts at that point are arbitrary since the 
Hamiltonian is independent of them. 

By applying these results to Eq. (20), it now follows that the structure 
of the optimum trajectories near cruise is given by the family of trajectories 
in Fig. (4). Specifically, no cruise segment occurs except at optimum cruise 

•wtgy * cope - 

By combining results from this and the preceding section, the Important 
result follows that, for the assumed engine model, the optimum trajectories, 
the corresponding optimum controls, and the performance are not affected 
by constraining the thrust to extreme values in the climb and descent segments. 


Case (B) ; Thrust Dependent 

A complete investigation of the neighborhood of the cruise point analogous 
to Case (a) requires estimates of the various thrust dependent derivatives in 
Eq. (26). However, understanding of this case can be obtained by examining 
the functions in Eq. (26) only along the thrust direction, i.e,, for AV ■ 0. 
Under that assumption, Eq. (26) simplifies to the following expression: 



<WS FC /V c ) ±l + (T c S FC /S FC ) + ( i AT l /2S 


FC ) I 2S FC + T c S FC , 

\ T t 2 / 


where the plus sign and AT > 0 are chosen for and the negative sign 

and AT < 0 for K. . 

an 
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This simplified approach focuses attention on the derivatives S and 

fc t 

S_„ , which are crucial for this case. The characteristics of these 

*C t 2 

derivatives can be deduced from plots of Sp C vs thrust* in Fig. 9. These 

plots* and those in Fig. 7* were derived from the operating instructions 

manual of a typical inservice turbofan (Ref. 6). Obviously* the assumption 

of a thrust-independent S pc is grossly violated for this engine* since at 

low thrust values the curves approach infinity; that is* they become 

undefined. However* at typical climb or cruise thrusts* corresponding to 

the upper half of the thrust range* the variation in S is only about 5 %. 

Fuel flow is also plotted in Fig. 9. The dashed line through the origin 

gives the best constant approximation to the fuel flow function. 

Comparison indicates an excellent match at high thrust* but an error of as 

much as 1200 lb/hr (550 kg/hr) at low thrust. For some applications the 

assumption of a constant Sp^ could be adequate if fuel flow errors at very 

low or idle thrust settings can be tolerated. 

For the upper two thirds of the thrust range* quadratic functions provide 

good fits to the Sp^ curves. Therefore, one can use the second-order 

Taylor series expansion at the cruise point to estimate Sp C for fairly 

large deviations of thrust from cruise thrust. 

The thrust in climb or cruise is typically larger than the thrust at 

which S__ is a minimum in Fig. 9. Both S_. and S will* therefore, 

rC ru.j, 

be greater than zero and so will the coefficient of AT in Eq. (32). It 

follows that the slope of Eq. (32) as a function of At is greater than zero 

for K and less than zero for K, . In other words, along the thrust 
up dn 

direction these functions have a strong minimum at the cruise point, whereas 
in Case (a) they were level along this direction. Along other directions. 
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the Investigation of Case (A) has shown a positive slope. Thus, if thrust 

is an unconstrained control variable along with airspeed, so that the cruise 

point lies in the Interior of the control region, then the optimum climb and 

descent thrusts and airspeeds will converge toward the optimum cruise thurst 

and airspeed as the climb and descent energies approach the cruise energy. 

It should be noted that his holds for all cruise energies, including those 

less than the optimum cruise energy, E . . Since the Hamiltonian is again 

cope 

zero at the cruise energy, it follows that the structure of the optimum 
trajectories as a function of range is identical to that of Case (A) ami is 
Illustrated by Fig. 4. Computer calculations for this case in Ref. 2, using 
a similar engine model, showed that the thrust is either maximum or idle 
for about three-fourths of the energy range between initial and cruise 
energies and then departs from the extremum values so as to converge smoothly 
to the value at cruise as cruise energy is approached. 

Consider now the case where thrust is constrained to some maximum in 
climb and is idle in descent. In that case, the minimum at the cruise point 
is not accessible since it does not lie in the region of permissible controls. 
Also, unlike Case (A), the thrust dependence of K and K^ n in Eq. (23) 
does not disappear along the thrust direction at V ■ V . Therefore, it is 
unlikely that at the minimum the sum of the two terms will be zero. The 
Hamiltonian is, in fact, greater than zero at any cruise energy. In order 
to show this, note in Fig. 9 that as thrust decreases, S FC increases without 
bound. It follows that 1^ will be less negative than it would be if S FC 
were thrust independent and therefore will be insufficient to cancel I at 
cruise energy, resulting in a positive value for the Hamiltonian. This was 
shown earlier to give rise to nonzero cruise segments below the optimum 
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cruls« energy. Thus, the structure of the optimum trajectories for the 
constrained thrust case Is given by the family of trajectories In Fig. 3. 

VII. Some Numerical Results 

Calculations of various types of optimum trajectories (minimum fuel and 
cost) have been carried out for the Boeing 727*200 and the McDonnell Douglas 
DC-10 aircraft using a digital computer implementation of the algorithm 
described in this paper. The models used in this program contain detailed 
propulsion and aerodynamic data, including thrust dependent specific fuel 
consumption and separate idle thrust and fuel flow models. The program also 
compensates for weight loss due to fuel burn as well as allowing the user to 
choose between the constrained and unconstrained thrust cases. The program 
is designed to augment or be a substitute for the computerized flight 
planning systems now used by airlines. A complete discussion of the results 
obtained in these calculations is beyond the scope of this paper. As an 
example of the results, Table 1 gives summary data for two minimum fuel and 
two minimum DOC trajectories for the 727-200. 

For the case of minimum fuel performance, the difference in fuel 
consumption between the constrained and free thrust optimum trajectories is 
63 lb (29 kg), or about 1%. For the case of minimum DOC performance, the 
difference is $0.03/n. mi., or again about IX. In both cases the differences 
have been found consistently, though dependent on range and aircraft weight. 
These differences, while not large, are about at the threshold (IX level) 
where they are considered significant in airline operating economics 

For both performance functions, the optimum cruise distances for the 
constrained thrust cases are a small percentage (11 and 182) of the total 
range. These results indicate that, for this engine, the specific fuel 
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consumption dependence on thrust does not have a strong influence on the 
trajectories. The IKK] penalty of flying minimum fuel trajectories is seen 
to be 14c/n. mi., or 6%, while the fuel penalty for flying minimum DOC 
trajectories is about 300 lb (136 kg), or 7Z. Each airline must evaluate 
the significance of the fuel and DOC differences between the trajectories 
in light of its schedules and route structures. 

The fuel cost and time cost factors used in these calculations are 
estimates obtained from a United States airline in mid-1977. These factors 
vary with time and between airlines. 

A comparison can also be made between the performance of optimum and 
currently used procedures. Such a comparison for a model of the DC-10 
flying a 220 n. mi. range yielded a 7 X fuel savings for the minimum fuel 
trajectories (Ref. 7). 


VIII. CONCLUSIONS 

The approach to trajectory optimization presented here has led to a 
rather detailed understanding of the characteristics of the optimum 
trajectories. The approach also lends itself to a numerically stable computer 
Implementation that can be incorporated in an airline flight planning system 
or, ideally, in an onboard performance management system. Furthermore, 
trajectories generated by this method can serve as benchmarks for evaluating 
other (suboptimum) algorithms. This possibility is especially intriguing at 
this time in view of the strong current effort in industry to develop the 
so-called On Board Optimum Performance Computers. 

Two pairs of opposing assumptions, constrained vs free thrust and 
dependence vs independence of specific fuel consumption on thrugt, played 
pivotal roles in determining the characteristics of the optimum trajectories. 
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If Che assumption of specific fuel consumption independent of thrust is 
justified, constrained thrust trajectories are identical to free thrust 
trajectories in structure and performance for the minimum fuel case. 

However, when the realistic dependence of specific fuel consumption on thrust 
is taken into account, there will be a difference in both performance and 
structure between constrained and free thrust cases. Tte actual numerical 
differences in performance depend on the propulsion and aerodynamic models 
as veil as otter factors and must be determined by computer calculations. 


Appendix 

It is to be proved that the loci of - AV ■ 0 and T - D ■ 0 are 

tangent at the cruise point, assuming that the cruise point at T * T , 

V ■ V is a minimum of the cruise cost W,/V along the locus T - D - 0. 

• c t 

This is equivalent to proving that the cruise point lies on both loci and 
that the slopes of the loci are identical at that point. 

That the cruise point satisfies - AV ■ 0 follows from the sequence 
of relations below: 



To prove that the slopes are identical, compute the gradient of 


W, 


AV 



- AV) 


A 

-IT 




+ j [ TS FC t + S Fcl 
l T J T -T 

f 

V-V 


(Al) 
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The perpendicular unit vectors i and j point in the speed and thrust 
directions, respectively. Now write X as a function of the perturbation 
AV: 


X » [(T +D AVJS^CT +D AV,V +AV)]/(V + AV) (A2) 

c v FC c v e c 


Since, by assumption, X has a minimum at V » V^, set the derivative of X 
with respect to AV equal to aero. This yields the following relation: 

X - VFC + T c( S FC t D v + S FC V ) ■ T c S FC /V c (A3) 

Next, compute the gradient of (T - D)(V/W) at the cruise point 


V(T - D)(V/W^ t - (V c /W)[i(-D v ) + j] 


<A4) 


v-v 


The slope of (Al) relative to the i direction is given by 
Siope = (t c 5 FCt + s,^/jr c S FCv - <T c s FC /v c )J 


(A5) 


After substituting Eq. (A3) in place of T c S p C ^ v in E< 1* (A5)» the slope 
simplifies to -1/D^, which is identical to the slope of Eq. (A4). 
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Table 1 Optimum trajectories for 200 n. mi. range, 

Boeing 727-200, JT8D-7 engines 
1965 United States standard atmosphere 

Initial weight: 136,000 lb (61,750 kg). Initial and final altitudes: sea level each, 

initial and final airspeeds: 210 knots each. 


Cost 

function 

Thrust constraint 

Cruise or 
maximum 
altitude 
(ft) 

Cruise 
distance 
(n. mi.) 

Fuel 

lib (kg)] 

Time 

(min:$) 

DOC a 

($/n. mi.) 

Minimum 

fuel 

Climb power in 
ascent, idle in 
descent 

29,500 

34 

4570 (2072) 

37:11 

2.54 


No constraints except 
operating limits 

31,300 

0 

4507 (2044) 

36:40 

2.51 

Minimum 

DOC 3 

Climb power in 
ascent, idle in 
descent 

30,400 

22 

4909 (2226) 

30:01 

2.40 


No constraints except 
operating limits 

30,500 

0 

4812 (2183) 

29:47 

2.37 

a Based on 

Cf - $0. 056/lb of fuel 

and ■ 

$410/hr. 
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Fig. 1 Assumed structure of optimum trajectories. 
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Fig. 3 Energy \/s range, H > 0 at E^. 


* 


31 



RANGE 


Kig. U Energy vs range, H - 0 at . 
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Fig. 6 Loci of T - D and W 
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Fir. H Dependence of operand function on dR and 


g at cruise enerfty 


E . 
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